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The electronic band structures of Bii_vSbx thin films can be varied as a function of temperature, pressure, 
stoichiometry, film thickness and growth orientation. We here show how different anisotropic single-Dirac- 
cones can be constructed in a Bii_fSbj^ thin film for different applications or research purposes. For predicting 
anisotropic single-Dirac-cones, we have developed an iterative-two-dimensional-two-band model to get a con- 
sistent inverse-effective-mass-tensor and band-gap, which can be used in a general two-dimensional system that 
has a non-parabolic dispersion relation as in a Bii_ jSbv thin film system. 
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Materials with two-dimensional (2D) Dirac cones in their 
electronic band structures have recently attracted considerable 
attention. Many important novel physical studies have been 
carried out on both massless and massive 2D Dirac fermions, 
including studies of the half-integer quantum Hall effect iH, 
the anomalous absence of back scattering, the Klein paradox 
effect iQl, high temperature superconductivity [jst] and unusual 
microwave response effects j^. The ultrahigh carrier mobil- 
ity of the Dirac fermions in graphene offers new opportuni- 
ties for a variety of electronics applications Isfc. Recently, 2D 
Dirac cones observed in topological insulators have identified 
this class of materials as candidates for quantum computation, 
spintronics, novel superconductors and thermoelectrics 

Materials with 2D single-Dirac-cones, especially 2D 
anisotropic single-Dirac-cones, are of special interest. Sim- 
ulations with ultracold atoms trapped on optical lattices 
have been used to study general 2D single-Dirac-cones |01. 
Graphene superlattice materials, are described by anisotropic 
2D single-Dirac-cones, which could potentially be developed 
for use in nano-electronic-circuits without cutting processes, 
and in table-top experiments that simulate high-energy rela- 
tivistic particles propagating in anisotropic space lUt]. 

Bii_,Sbv alloys could be especially attractive for studies 
of 2D Dirac fermions. Firstly, the band structure of three- 
dimensional (3D) bulk Bii_vSbv can be varied as a function 
of Sb composition x, temperature, pressure and strain, and the 
Fermi level can be adjusted to change the electronic properties 
Jgt]. Secondly, not only have bulk state Dirac points been stud- 
ied in bulk Bii_xSbv alloys ifioll . but also the first observation 
of a surface state 2D single-Dirac-cone for a topological in- 
sulator was made in bulk Bio.gSbo.i 111 ill . Some experiments 
on Bii_vSbv thin films grown normal to the trigonal direc- 
tion have already been carried out Il2l - ll4ll . However, possible 
Dirac-cone materials of 2D Bii_vSb;c thin films have not yet 
been studied. 

In this Letter, we have systematically studied how the 
properties of the 2D L-point Dirac cones in Bii_rSb, thin 
films, such as their anisotropy and linearity, change with film 
thickness and growth orientation. We have also discovered 



that different varieties of Dirac cones can be constructed in 
Bii_xSbv thin films, such as single-Dirac-cones with differ- 
ent anisotropies. For the methodology, we have developed an 
iterative-two-dimensional-two-band model, which could be 
used as a general model to study systems with non-parabolic 
2D dispersion relations or Dirac cones. 

We recall that the band-gap and band-overlap of bulk Bi 
can be changed by adding Sb to form Bii_,Sbv alloys. 3D 
Dirac points may form under proper conditions at the three 
L points ifioll . where the valence band and conduction band 
exhibit a band-crossing. At each L point, ^(k) becomes lin- 
ear in the limit of band-crossing {Eg — ^ 0), since ^(k) ~ 

± ( (v ■ p)^ + E^) - , where v is the carrier velocity, p is the car - 



rier momentum, and Eg is the L-point direct band-gap 11 1LI15II . 
Historically, the 3D two-band model |fl6ll has been successful 
in describing the non-paraboHc dispersion relation £ (k) in 3D 
bulk Bi materials by 



p-a-p = £(k)(l + ^), 



(1) 



where a is the inverse-effective-mass-tensor, and we assume 
here that a is the same for both the conduction band and va- 
lence band within the context of a two-band model and strong 
interband coupling. Generally, the relation between a and E^ 
around an L point is described as 
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under the k ■ p approximation IitIi . where I is the identity ma- 
trix and mo is the free electron mass. We assume here that 
a two-band model also applies to Bii_^Sb^ alloys, where the 
influence on a of adding Sb atoms up to a Sb concentration of 
x=Q. 1 to bulk Bi iH is then 
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We developed an iterative-two-dimensional-two-band 
model for Bii-^Sb, thin films, where both a^'''" and Eg''"' of 
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a Bii_.,Sb.v thin film differ from a*"'*^ and E^""' in the bulk 
case, and they are both unknown at the start. The conduction 
band and valence band at an L point are firstly separated 
by sf^ = E^^""^{Bi) and the inverse-effective-mass-tensor is 
„[0] ^ a'"'"'{Bi). In a Bii_;,Sb, thin film, the L-point gap will 
be influenced by both the added Sb atoms and the quantum 
confinement effect, so that the L-point gap to lowest order 
becomes 
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The change in the L-point gap will lead to a change of the 
L-point inverse-effective-mass-tensor according to 
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where n denotes the step number in the iteration. Eq. (5) is a 
general equation for iterations, which is just a consequence of 
Eq. (3). Thereafter, the L-point band-gap can be updated with 
the new inverse-effective-mass-tensor by 
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The procedure is repeated until a'"' and E™ become self- 
consistent, and then we get an accurate solutions for al"' = 
a^''"'{Bii_^Sb^) and E^ = E{"'"{Bii^^Sb^) for the Bii^.^Sb,- 
thin film. The Hamiltonian for Bi and Bii_^Sbv based on k • p 
theory in this model is equivalent to a Dirac Hamiltonian with 
a scaled canonical conjugate momentum 11511 . Thus, Eqs. (5) 
and (6) are also good approximations to describe the Dirac 
cones. The band parameters we use in the present calcula- 
tions are values that were measured by cyclotron resonance 
experiments 1IT9I1 . 

In order to get a single-Dirac-cone, where only one Dirac 
cone effectively contributes Dirac carriers in the 2D system, 
we need to grow the Bii-vShx thin film normal to a low sym- 
metry direction so that we can have a single L-point that dif- 
fers from the other two by breaking the 3 -fold symmetry oc- 
curring in 3D. Moreover, two quantities in Eq. (4) need to be 
minimized, namely the bulk-band-gap term E'^"'^{Bi\_^Sbx) 

and the quantum-confinement-induced term h^u}^\/A-l} for 
this special single L-point. Thus, by growing a film normal 
to the bisectrix direction the 3-fold symmetry is broken and 
the inverse mass component is near its minimum lfl9ll . 
which is discussed below. Minimizing the value of bulk term 
£''"''*(B;i_v5'i'.v) by varying temperature T, pressure P and Sb 



composition X has been discussed in literature III iLIITllISll . es- 
pecially for the case of T=77 K P=l atm and x=0.04, where 
£'^"'*'(B;i_.v5/7.v) ceases to 0. Bii_;rSb,^ thin films with this 
composition of x=0.04 (Bio.96Sbo.04) have been already syn- 
thesized in experiments by researchers in Ukraine jioll for 
films grown normal to the trigonal direction. In the present 
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FIG. 1: (color online). Schematics of the conduction band and va- 
lence band for different carrier pockets of a Bio.96Sbo.04 thin film vs. 
film thickness: (a) illustrates the T-point hole pocket and the three 
L-point electron pockets in the first Brillouin zone for bulk bismuth 
with a rhombohedral crystal structure, (b)-(d) are for Bio,96Sbo.o4 
thin films that are grown normal to the bisectrix direction associated 
with the lC point. E=0 corresponds to the middle point of the L- 
point band-gap. (b) and (c) are comparisons between energies of the 
conduction band at the three L points, and l'"'' for different 

film thicknesses. AL^^''^'^' is smaller than the thermal smearing the 
energy of a carrier at at 77 K which corresponds to < 310nm 
from the iteration of Eqs. (4), (5) and (6). (d) Comparison between 
energies of the conduction band and valence band at the L''^ point 
and the valence band at the T point, (e) and (f) are comparison be- 
tween the energies of the conduction band at the three L points, L^'', 
L(-) and l'-'', for different film thicknesses /- when the Bio,96Sbo,o4 
thin film is grown normal to the [6061] direction. 



paper, we illustrate simple examples of the various kinds of 
single-Dirac-cones that can be found in Bio.96Sbo.04 for more 
generally and in particular for conditions where Dirac cones 
can be found. 

According to our calculations, a single 2D Dirac cone can 
be constructed in a Bio.96Sbo.04 thin film by breaking the 3- 
fold symmetry of the L points. We calculate here the lowest 
conduction band and the highest valence band at the three L 
points, namely l''', L*^^' and l'^' shown in Fig. 1(a), when 
the Bio.96Sbo.o4 thin film is grown normal to the bisectrix 
direction associated with the L''' point [Fig. 1(b)] in this 
case. The L'^' and L'^' points still retain a double degen- 
eracy with each other, while the energy difference between 
the lowest conduction band at the L''' point and the corre- 
sponding energy at the L*^^' point or L'^' point (AL^^''^'^' = 
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) increases as the film thick- 



ness decreases [Fig. 1(c)]. Below 77 K the band structure of 
Bii_;cSb^ does not vary significantly with temperature. At 77 
K, the smearing of the Fermi-Dirac distribution is ^ A:bL=6.64 
meV. Thus, when AlI^''^'"' > 6.64 meV, the Dirac fermions 



can all be associated with the L'^'^ point. The energy differ- 
ence A^l''^'^' = 6.64 meV corresponds to a film thickness 
of h = 310 nm as shown in Fig. 1(c). Thus, we can have a 2D 
single-Dirac-cone at the L^'^ point in a Bio.96Sbo.04 film that is 
grown normal to the bisectrix direction and thinner than 310 
nm. 

In experiments, a Dirac cone is usually not perfect. A mini- 
gap often exists which induces a tiny mass at the apex of 
the Dirac cone. Such an effect also occurs in single layer 
graphene There are two main features that characterize 
the properties of a Dirac cone, the tiny mass occuring at the 
apex of the Dirac cone and the fermion velocity v(k) as a func- 
tion of k. Because the fermions are Unearly dispersed near a 
Dirac cone, v(k) should only be a function of the direction of 
(k). The anisotropy of the Dirac cone can be characterized by 
the ratio between the maximum and minimum v(k). 

The Dirac cone at the L^'' point in a 310 nm thick Bii_YSb;c 
thin film grown normal to the bisectrix axis has a linear ^(k) 
behavior with a very sharp apex [Fig. 2(a)] . In this film, k, and 
ky are in the trigonal direction and the binary direction, respec- 
tively. Ej}^ for this Dirac cone is smaller than 0.1 meV, and 
the effective mass at the apex of the Dirac cone is ~ 10^^;7Jo, 
which can be considered as essentially massless. We have also 
calculated the v(k) relation of the Dirac fermions for different 
values of momentum [Fig. 2(b)]. It can be seen that the max- 
imum and minimum of v(k) are 1.6 •co/300 (along kv) and 
1.1 ■ Co/300, (along ky) respectively, which differs by a factor 
of ~^1.5, where cq is the speed of light. 

In some applications, we need all the carriers within the 
film to be linearly dispersed, which also means that the Fermi 
level Ef needs to be arranged to have no T-point paraboli- 
cally dispersed carriers. This condition is not met in bulk 
Bio.96Sbo.o4 , but it can be achieved in a Bio.96Sbo.04 thin film. 
The highest valence band at the T point by comparison to that 
of the lowest conduction band at the l''' point is calculated 
in Fig. 1(d). Based on Fig. 1(d), when /, < 40 nm the film 
becomes an indirect gap semiconductor and if is adjusted 
to be high enough above the top of the valence band, carriers 
can all be L'^'^-point electrons, with no T-point holes. At this 
film thickness, Eg^ '' at the l''^ point is smaller than 4 meV, and 
the effective mass at the apex of the Dirac cone is ~ lO^^wzo, 
which is still essentially massless. The corresponding Dirac 
cone is plotted [Fig. 2(c)], as well as the velocity vs. mo- 
mentum relation v(k) for the Dirac fermions [Fig. 2(d)]. The 
linearity and the anisotropy of the Dirac cone is not notably 
influenced by film thickness for a film of /- = 40 nm, as can 
be seen by comparing Figs. 2(c) and 2(d) to Figs. 2(a) and 
2(b), respectively. When h<3l nm, the film become a direct 
gap semiconductor [Fig. 1(d)], and the both the Dirac elec- 
trons and holes can be free from the influence of the T-point 
parabolically dispersed holes if the Fermi level is above the 
top of the valence band. 

In some other applications, a higher anisotropy of the Dirac 
cone could be required, e.g. in nano-electronic-circuit design, 
^(k) at the L'^' point of a Bio.96Sbo.04 thin film grown nor- 
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FIG. 2: (color online). Different anisotropic single-Dirac-cones in 
different Bio.ci6Sbo.o4 thin films for various applications or research 
purposes: (a) and (b) describe a sharp-apex anisotropic single-Dirac- 
cone in a 310 nm thick Bio.96Sbo.04 film grown normal to the bisec- 
trix direction. For convenience, the origin of momentum k is chosen 
to be at the point, (c) and (d) describe an anisotropic single- 
Dirac-cone where the Dirac electrons are free from the influence 
of r-point holes in a 40 nm thick Bio.96Sbo.04 film grown normal 
to the bisectrix direction, (e) and (f) describe a highly anisotropic 
single-Dirac-cone in a 300 nm thick Bio.96Sbo.04 film grown normal 
to the [6061] direction, (a), (c) and (e) show the dispersion relation 
of single-Dirac-cones. (b), (d) and (f) show the group velocity v of 
Dirac fermions at different values of momentum k. (c)-(d) are not 
significantly different from (a)-(b), and (e)-(f) are obviously more 
anisotropic than (a)-(b) and (c)-(d). 



mal to the corresponding bisectrix direction has a smaller tiny- 
apex-mass but a lower anisotropy, while ^(k) at the L''^ point 
of a Bio.96Sbo.o4 thin film grown normal to the trigonal di- 
rection has a larger tiny-apex-mass but a higher anisotropy, 
when compared to Bio.96Sbo.04 thin films of the same thick- 
ness. To get a single-Dirac-cone with both a small tiny-apex- 
mass and high anisotropy, we can choose to grow the film 
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normal to a direction between the trigonal [0001] axis and the 
bisectrix axis [1010] within the binary plane, where [0001] and 
[lOlO] are hexagonal notations for such crystalline directions 
in bulk rhombohedral Bii_;fSb^. As an example, we illus- 
trate in Fig. 2(e) and (f) the film grown normal to the [6061] 
direction, which is a lower-symmetry direction in the binary 
plane with an angle of 14° to the bisectrix axis and 76° to the 
trigonal axis. Results for the lowest conduction band and the 
highest valence band at the three L points as calculated for 
a Bio.96Sbo.o4 thin film grown normal to the [6061] direction 
are shown in Fig. 1(e). According to the criterion we used 

above [Fig. 1(b)], the energy dependance between i?!'' and 
is A£|''^'^' < 6.64 meV when I- < 300 nm [Fig. 1(f)]. 
The l'') point Dirac cone in the 300 nm thick Bio.96Sbo.o4 
thin film grown normal to the [6061] direction is plotted [Fig. 
2(e)]. In this film, ky is still normal to the binary direction, 
and k.v is in the direction perpendicular to both the binary di- 
rection and the [6061] direction. Eg^^ for this Dirac cone [Fig. 
2(e)] is smaller than 0.46 meV, and the effective mass at the 
apex of the Dirac cone is still negligible (~ lO^^mo). The 
v(k) relation of the Dirac fermions calculated for this cone is 
shown in Fig. 2(f). The linearity and the sharpness of ^(k) 
are not sacrificed significantly compared to Fig. 2(a). How- 
ever, the anisotropy of the Dirac cone in Fig. 2(e) is increased 
remarkably. The maximum and minimum values of the Dirac 
fermion velocity v(k) near this Dirac cone are, respectively, 
1.65 ■ Co/300 and 0.55 • co/300, which differ by a factor of 
~3 [Fig. 2(f)], which we consider to be a highly anisotropic 
single-Dirac-cone. 

The Fermi level can be adjusted by doping such a Bii_vSbr 
sample with foreign elements, as for example by doping the 
sample with selenium or tellurium (Column VI) to increase 
the Fermi level (n-type) and by doping with tin or lead (Col- 
umn IV) to decrease the Fermi level (p-type). In the graphene 
system, the Fermi level cannot be adjusted freely by heavy 
doping, because the band distortion due to heavy doping may 
destroy the linear k aspect of the Dirac cone itself. However, 
in a Bii-vSbx thin film, the Sb composition, e.g. x=0.04, is 
over 2 orders of magnitude larger than normal electronic dop- 
ing concentrations (< 10^^), so that the band structures will 
not be obviously distorted by electronic doping. Thus, the 
Fermi level can be varied over a wide range of doping levels 
relative to £ = in a Bii_vSbv thin film. 

In conclusion, we have developed an iterative-two- 
dimensional-two-band model to describe a general 2D non- 
parabolic anisotropic dispersion relation. Based on this 
methodology, we have made a prediction that anisotropic 
single-Dirac-cones can be constructed in Bii_;cSb;i thin films. 
Some critical cases of L'^') point single-Dirac-cones are illus- 
trated as examples: a single-Dirac-cone with a sharp apex 
[Fig. 2(a)], a single-Dirac-cone without influence of T- 
point parabolically dispersed holes in a Bii_j.Sb;r thin film 
[Fig. 2(c)] and a single-Dirac-cone with high momentum 
anisotropy [Fig. 2(e)]. Novel physical phenomena associated 
with Dirac cones and Dirac fermions that have been previ- 



ously reported in other materials systems could hopefully also 
be observed in Bii_,Sbv thin films. Because the Bii_rSbv 
thin film system has special features as we discussed above, 
we can also expect to observe new physical phenomena that 
have never been observed in other systems. 
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